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In this paper we present a phase classification of (effectively) two-dimensional non-Abelian nemat- 
ics, obtained using the Hopf symmetry breaking formalism. In this formalism one exploits the 
underlying double symmetry which treats both ordinary and topological modes on equal footing, 
i.e. as representations of a single (non-Abelian) Hopf symmetry. The method introduced in the 
literaturei*- and further developed in a paper published in parallel- allows for a full classification of 
defect mediated as well as ordinary symmetry breaking patterns and a description of the resulting 
confinement and/or liberation phenomena. After a summary of the formalism, we determine the 
double symmetries for tetrahedral, octahedral and icosahedral nematics and their representations. 
Subsequently the breaking patterns which follow from the formation of admissible defect conden- 
sates are analyzed systematically. This leads to a host of new (quantum and classical) nematic 
phases. Our result consists of a listing of condensates, with the corresponding intermediate residual 
symmetry algebra % and the symmetry algebra U characterizing the effective "low energy" theory 
of surviving unconfined and liberated degrees of freedom in the broken phase. The results suggest 
that the formalism is applicable to a wide variety of two dimensional quantum fluids, crystals and 
liquid crystals. 
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I. INTRODUCTION 



The goal of this paper is to apply a Hopf symmetry 
breaking analysis to defect condensates in nematics. The 
subject of quantum liquid crystal phases is quite exten- 
sive, and we first highlight some nematic phases where 
our methods are the most relevant. Then we describe 
how Hopf symmetries (in particular, double symmetries) 
characterize the degrees of freedom in phases with spon- 
taneously broken symmetry, and discuss the Hopf sym- 
metries relevant to the exotic nematic phases we focus on 
in this work. Finally, we describe the formalism for sym- 
metry breaking, and work out all possible defect medi- 
ated phase transitions in exotic achiral nematics, be they 
classical or quantum. Such phases have only been ob- 
served in "classical" systems, in which we do not expect 
quantum superpositions of defects to form condensates. 
Still, we choose these phases because they highlight the 
power and generality of our approach. There are known 
examples of quantum liquid crystals in which a rank two 
tensor order parameter field is needed to describe the 
phase, and our methods are definitely applicable to such 
systems, though capturing the exotic nonabelian phases 
requires at least a third rank tensor. 
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A. Defect condensates 

Classical liquid crystals have been studied for a long 
time. Recently there has been a renewed interest in ex- 
otic liquid crystals, and they have been invoked to explain 
certain phases in bent-core liquid crystals. An exhaustive 
analysis of ordinary symmetry breaking patterns in clas- 
sical nematics can be found in the literature^. In a sense, 
our study complements this analysis, as we work out an 
exhaustive analysis of defect mediated phase transitions 
and their interpretation as the breaking of certain double 
symmetries. 

The literature on liquid crystal phases in quantum Hall 
systems has vastly grown in recent years^. In High-Tc 
superconductivity, the stripe phase is a two dimensional 
quantum smectic, and recently a theoretical study has 
analyzed the possibility of having a topological nematic 
phase in such a system^. The nematic order is arrived at 
by a defect condensate from a crystalline phase, in other 
words through a non trivial vacuum expectation value 
of some disorder parameter. This in contrast with the 
more conventional spin nematic order, known to exist 
in superfluids* and High-Tc superconductors^, where a 
biaxial nematic phase has been found. Research on spin 
nematic order is still vibrant. 

Our work offers a general approach to the study of all 
conceivable condensate mediated phase transitions using 
an analysis similar in vein to Landau theory. Just like 
Landau theory, it could serve as the kinematical back- 
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bone of more detailed dynamical studies (involving effec- 
tive Hamiltonians and a renormalization group analysis, 
for example). In this paper we focus on exotic nemat- 
ics whose residual symmetry correspond to the tetrahe- 
dral, octahedral or icosahedral group. We choose such 
phases to expose the power of our method, and because 
as described above there is an ever growing interest in 
quantum liquid crystal phases. 

II. NEMATIC PHASES AND DOUBLE 
SYMMETRIES 

In this section we define nematic phases quite gener- 
ally in terms of symmetries, focusing on systems that are 
effectively two dimensional. We then briefly recall how 
the representation theory of an underlying Hopf (quan- 
tum double) symmetry leads to a description that treats 
regular excitations, topological point defects and dyons 
on equal footing (for details, see a paper published in 
parallel"^). On rather general grounds we determine the 
Hopf symmetries that characterize the relevant nematic 
phases. The precise outcome for nematics is a modified 
quantum double, which is a variation on Drinfeld's quan- 
tum double of a groupSS. The generality of the approach 
makes our methods applicable to basically all "nematics" , 
be they classical or quantum, global or gauged. In fact, 
any phase that results from spontaneous condensation 
phenomena can be subjected to a similar analysis. Our 
first task is to classify all modes in these media as irreps 
of a Hopf symmetry. This classification serves as a cru- 
cial ingredient in the analysis of the next section, where 
we study phase transitions induced by the condensation 
of defects. 



A. Nematic liquid crystals 

A nematic liquid crystal is a phase with com- 
plete translational symmetry, and incomplete rotational 
symmetry^*^'^^'^^. The phase then inherits the name of 
the residual rotational group: if the residual rotational 
group is the tetrahedral group, for example, the phase is 
called a tetrahedral nematic. 

The residual internal rotational symmetry group H can 
be any proper subgroup of G = 5*0(3). H is the stabi- 
lizer of some fixed tensor in a representation of G. This 
implies that H must be a closed subgroup of 5*0(3). The 
closed subgroups of 50(3) are well known: 

G = 50(3) ^He {Gn, Dn, T, O, /, 50(2) >^ Z2}, (1) 

where we used the notation employed in the crystallog- 
raphy literaturei^. 0„ = Z„, the abelian cyclic group 
of order n. Dn is the dihedral group of order n, T is 
the tetrahedral group, O the octahedral group, and / the 
icosahedral group. 

If we want to consider inversion symmetries, then we 
break G = 0(3) to a closed subgroup H of 0(3). In 



the case of an achiral tetrahedral nematic, 0(3) is bro- 
ken to Td (we adopt the crystallographic notation for the 
symmetries of a tetrahedron including reflections). 

We will call the residual symmetry group H the electric 
group of the phase, and in general we will denote an 
electric group of a phase by Hei ■ 

B. Excitations 

In general, in a phase where a group G is spontaneously 
broken to a subgroup H, one distinguishes between three 
types of modes: regular excitations in what is often called 
the "electric" sector , topological defects corresponding 
to the magnetic sector, and mixed excitations in the dy- 
onic sector. 



1. Regular excitations 

Regular excitations (or regular modes) are smooth, low 
energy excitations of the basic fields that characterize 
the system. Examples are continuity modes (present be- 
cause of conservation laws), and Goldstone modes. These 
modes may be coupled to each other (such as is the case 
in classical nematicsi^ ). We will often refer to these 
regular modes as electric modes. 

The regular modes transform under irreps H^ of the 
symmetry group H. The corresponding states form a vec- 
tor space on which the elements of H act as linear trans- 
formations. The states are denoted as >, where the 
^■^ stand for all the numbers we need to characterize the 
state, spatial coordinates and other quantum numbers. 
Multi-particle states are described by tensor products of 
the elementary representations which are assumed to be 
reducible and can be decomposed into irreducible com- 
ponents given by the appropriate Clebsch-Gordan series: 

n„ ® H/3 = iVj^H^ , (2) 

These rules for combining representations are also called 
fusion rules. They imply that if we bring particles 1 
and 2, in states > and |^2 > respectively, closely 
together, and we measure the quantum numbers of the 
combined system, that we can get different outcomes, 
precisely given by the Clebsch-Gordan coefficients corre- 
sponding to the decomposition Q. 

Let us consider one of the irreps H^ , and choose a basis 
for the corresponding vector space denoted by {|e" >}, 
where j labels the different basis vectors. We then write 
for the (matrix) action oi g € H in that basis: 

aig) ■ \e'^ a{g)'; • > . (3) 

The physical requirements on the representations are 
that they are unitary, because the state should stay nor- 
malized under the action of G. However, if H is not sim- 
ply connected, then in a quantum nematic the states can 
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transform under projective representations of H. This 
means that the action of H on the states is not a group 
homomorphism: the action of 92 first and then gi is not 
equal to the action of (7152 • The actions may only differ 
by a phase: 

a(5i)a(52)«' = e'^<^aug2)^(^g^g^^^, (4) 

This is allowed because the phase factor disappears when 
we calculate (transition) probabilities For exam- 

ple, half integral spin representations transform under a 
projective representation of SO{3). As a matter of fact 
projective representations of a group H correspond to 
(faithful) representations of H, defined as the lift of H 
in the universal covering group G of G. For example, a 
spin i particle forms a doublet, which is a faithful rep- 
resentation of SU{2). 

2. Topological defects 

Topological point defects in two spatial dimensions (or 
line defects in three dimensions) correspond to nontriv- 
ial configurations of some order parameter field which 
are stable for topological reasons. The point defects are 
characterized by the first homotopy group 'Ki{G/H) of 
the vacuum manifold G / H . The group element that cor- 
responds to a given defect is called its topological charge 
or magnetic flux. In general we will denote the magnetic 
group of a phase by Hm ■ 

Using a standard theorem from homotopy theory: 

7ri{G/H)^7roiH), (5) 

we find that the point defects are characterized by the ze- 
roth homotopy group (which studies the connected com- 
ponents) of H. In particular, if _ff is a discrete group, 
then Ili{G/H) = H. Let us assume that H is discrete, 
then we will call Tii[G/ H) the magnetic group Hm of the 
phase. The first observation we should make about the 
composition rule for the defect charges is that it is spec- 
ified by the structure of the first homotopy group and 
corresponds therefore to group multiplication. We have 
to say more about this though, because in the cases of 
interest these groups are non-Abelian which at first sight 
gives rise to unwanted ambiguities in the fusion rules for 
defects. 

Let us now denote the (internal) physical state corre- 
sponding to a defect charge 5 by a ket \g >. The defect 
states form a vector space V spanned by the \g > with 
c, e V = {Y,. \j\gj >: Xj £ C,gj £ H]^. The vec- 
tor space spanned by group elements and equipped with 
the group multiplication is called the group algebra and 
denoted CH. 

In our quantum treatment it can make sense to add 
certain defect states. A state l^i > +\g2 >, for exam- 
ple, would correspond to a quantum superposition of the 
defects gi and g2- A priori there is no obvious classical 
interpretation for this superposition. We note however 



that there are actually cases where superpositions of de- 
fects can be given a classical interpretationi^. 

States in different irreps cannot form a superposition: 
the span of states in one irrep forms a superselection 
sector. To figure out how many defects are in one sector, 
we need to know how the electric group H^i acts on the 
magnetic group Hm- As a matter of fact, if we have a 
defect in our system and act on the system with a global 
symmetry transformation, then we may obtain a different 
defect. Given h £ H^i and g £ Hm, we denote the action 
of h on g hy h ■ g. This action satisfies 

hi ■ (ft.2 • 5) = {hih2) ■ g V/ii, ft.2 £ H^ug £ H^ 
h ■ {9192) = [h ■ gi){h ■ 32) V/i £ Hei,gi,g2 £ Hm 

The first equation is natural, it simply says that the 
group Hei acts on Hm as a group. The second equation 
implies that the action of a global symmetry transforma- 
tion ft, on a configuration that is composed of two defects 
next to each other, with topological charges 91 and 32 , is 
equal to the action on each defect separately with h. 

For example, if Hm = Hei = H, the action of a global 
symmetry transformation h £ H on defect > is a group 
conjugation of the topological charge: 

h-\g>^ \hgh-' > (6) 

The defect representations are therefore labelled by the 
defect classes A in Hjyr^. These classes correspond to 
sets of defects that are transformed into each other un- 
der the action of H^i- We should think of the classes 
A as defect representations 11"*, and they represent H^i 
invariant sectors of the theory. 

At this stage of the analysis the defect representa- 
tions 11"* might be reducible. However, the algebra can 
be extended with other operators in our theory which 
make the classes irreducible representations. We can 
in principle measure the precise fiux of a defect, us- 
ing (global) Aharonov-Bohm scattering experimentsiSiii. 
Correspondingly there exist certain projection operators 
Pg in our theory [g £ Hm), that act on the defects ac- 
cording to 

Pg\g' >^ 5g,g,\g' > . (7) 

These projection operators span a vector space which 
is isomorphic to the vector space of functions from Hm 
to C, which we denote by F{Hm)- Namely Pg can be 
associated with the function on Hm defined by Pg{g') = 
^g.g'- F{Hm) Can be turned into an algebra by taking 
pointwise multiplication, and this is precisely the algebra 
of the projection operators! Thus we will associate the 
projection operators with F{Hm)- 

The defect classes form irreps under the combined ac- 
tion of Hei and F(Hm)- Note that we've described the 
action oiF{Hm) and H^i separately, and we need to know 
what happens when a projection operator and a global 
symmetry transformation are applied in succession. Thus 
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we want to turn the combination of H^i and F(Hm) into 
an algebra, i.e. we want to be able to multiply elements 
of Hei and F{H„i). Physics dictates what the answer is: 
the multiplication in this algebra is set bjsi^ 

hPy = Pn.gh. (8) 

The physical motivation for this equation is as follows: if 
we measure a flux g with Pg, and then conjugate the de- 
fect with h, we have a flux h ■ g. This action is equivalent 
to first acting on the defect with h, and then measuring 
h ■ g with Ph.g. 

We call the algebra defined in this way a modified quan- 
tum double (because it closely resembles the quantum 
double D{H), which is a special case of the modified 
quantum double with H„i ~ Hei = H), and denote 
it bySi F{Hm) x CHei- As a vector space, the mod- 
ified quantum double is simply F{Hm) 'S) CHgi. The 
multiplication is set by the action defined above. Thus 
we conclude that the defects transform under irreps of 

F{Hra) X <CH,i. 

The tensor product \gi > ®\g2 > of two defects is to be 
interpreted as "a configuration with defect gi to the left 
of defect 172" ■ The order is important: if we measure the 
total flux of \gi > \g2 > we get gig2, while the total flux 
of \g2 > \gi > is 3231- Thus we define the action of the 
projection operators on the tensor product as follows: 

Ph{\gi > «)|52 >) = Sh,g,gA\9l > '»|52 >)• (9) 

We now give a couple of examples of nematic phases 
and their associated modified quantum double: 

• A chiral tetrahedral nematic 

A tetrahedral nematic has internal symmetry 
Hei = T, where T is the tetrahedral group. There 
are no reflections in H because the phase is chiral. 
The magnetic group is i?™ = ni(S'Cl(3)/r) = f, 
where T is the double cover of T in SU{2). 
Therefore the relevant modified quantum double 
is ^ = F{T) X T. If we are considering a quantum 
mechanical nematic, and we have spinors around, 
then A = F{f) x f ee D{f). This is the quantum 
double of f . 

• An achiral tetrahedral nematic 

Now Hei — Td, the group of symmetries of a 
tetrahedron including reflections. The magnetic 
group is the same as in the chiral case: H„i — 
ni(0(3)/Td) = ni{SO{3)/T) = f. Therefore 
A = F{T) X Td- If we allow for spinor electric 
irreps, thei^ A = F{f) x f^. 
Note that the algebra multiplication is determined 
by the action of the electric group on the magnetic 
group, which we still need to calculate. We do this 
in the next section. 

• A uniaxial nematic 

The local symmetry is Hei = S0{2) xi Z2. The 



magnetic groups is ni(S'C>(3)/S'0(2) XZ2) ~ Z, but 
the analysis is actually more subtle: the relevant 
modified quantum double turns out to be 

A = F(S'0'(2fx Z2) X CSO(2) x Z2. (10) 

where S0{2) x Z2 is the cover of 50(2) x Z2 in 
SU{2). The defects carry a continuous label. We 
stated before that Hm is always discrete, and in- 
deed in this case Ui{SO{3)/SO{2) x Z2) ~ Z2. 
Thus there is only one nontrivial defect homotopy 
class. Given a configuration of the fields that cor- 
responds to this defect, we can always act on the 
configuration with the residual rotational symme- 
tries of SO{2) C Hei to obtain a continuous family 
of configurations that correspond to the same de- 
fect. Because the defects in this family all have the 
same energy, we should expect a zero mode in this 
defect sector. This in turn leads to interesting phe- 
nomena, such as Cheshire charge (for a discussion, 
in a two dimensional context, see referencesii^) . 
Cheshire charge is usually associated with gauge 
theories, but it can exist in global theories as welUS. 
We point out that the Hopf symmetry approach is 
apphcable to these cases, though beyond the scope 
of this paper. 



3. Dyonic modes 

We now want to complete our description of the full 
"internal" Hilbert space by including the mixed sectors 
carrying both nontrivial electric and magnetic charges. 
For every defect class A we choose a preferred element 
gA as a representative, then all defects can be written 
a.s h ■ gA for some A and some h e H,n. Call Na the 
normalizer'^'^ of gA- The normalizers of elements in the 
same defect class are isomorphic: Nh-g^ = hNAh~^. 
A dyonic mode is an electric mode in a topologically non- 
trivial sector corresponding to a defect gA- In that case 
there is an important restriction due to the topological 
obstruction to globally implement all global symmetry 
transformations of Hei^^-^^. Only the subgroup Na can 
be globally implemented, and hence the electric modes 
in such a sector will transform under an irrep a of Na - 
Extending our ket notation to all dyonic/mixed sectors 
we denote a state with electric component |e" > and a 
defect h € A in the background by \h,e'^ > (following 
the notation in the literature^^). The |e" > form a basis 
of the vector space on which a acts, so that the e" > 
(with h G A) are a basis of the vector space associated 
to the dyon. We denote this irrep of our dyon by U^. 
The action of global symmetry transformations on this 
vector space is subtle. If we take a transformation 
n G Na, then 

n-\h,e°' >=niin)\h,ef >= \h,a{n)ef >, (11) 
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i.e. n acts on the electric mode. 

But if the transformation g ^ Nj^, it will transform the 
defect, while at the same time it can act on the electric 
mode! To describe this action, it is convenient to define 
another notation for the vectors in 11^. 
First note that the elements of the defect class A are 
in one-to-one correspondence with left iVyi-cosets in Hei ■ 
Choose representatives xf of left Na cosets, such that 
= e. Then xf corresponds to hf = xf ■ gA, where 
hf is an element of A. This association is well defined 
because it is independent of the particular choice of rep- 
resentative xf of the left Na coset, since by definition the 
elements of Na commute with gA- Furthermore, different 
xf correspond to different elements hf of A, and we have 
A = {hf = gA, hf, hf, . . . , hf} . Now a basis of the vec- 
tor space on which Hf acts is given by {\hf,e'j >}. 

Alternately, we can denote \hf ,6" > hy xf x le" >. In 
this notation, acting on the defect hf with g corresponds 
to multiplying xf by g from the left. Thus 

g-\hf,e'i^ >=gxfx\e'i^ >=xfnx\e'i^ >= \hf,a{n)e'^ > 

where gxf = xfn, with n e Na- In other words, gxf sits 
in some left Na coset. Since the xf form representatives 
of left Na cosets, gxf is equal to xfn for some k and 
some n £ Na- This n then acts on the electric part of 
the dyon. This notation is the most transparent notation 
we can adopt for the action of g on the dyon. 
The action of the projection operator Ph on the dyon is 

PH-\hf,e^>=6,,^^A\hf,ef>, (12) 

thus it projects the defect part. 

Summarizing, the 11^ are irreps of F{Hm) x CHei- It 
turns out that these are all the irreps^'* of F{Hm) x CHei- 
We denote the vectors on which 11;^ acts by V^- 

Note that the electric and magnetic modes discussed 
are also irreps of F{H„i) x CH^j- Namely, electric 
modes are irreps H^", with Ce the conjugacy class of the 
identity e: Ce — {e}, and the excitations carry irreps of 
the full group, i.e. Na = H^i- Magnetic modes are irreps 
(where 1 is the identity or trivial representation), 
dyons with a trivial representation of Na- In this 
sense the quantum double offers a unified description of 
electric, magnetic and dyonic modes- 

The steps towards classifying all the irreps of the dou- 
ble F{H„i) X CHei are relatively straightforward^: 

• Determine the defect classes A of iJ„i, i.e. the 
classes under the action of H^i 

• Pick a preferred element gA for every A, and deter- 
mine the normalizer Na of gA 

• Determine the irreps a of Na 

• The hrcps of F(iJ,„) x CH^i are the set {llf}- 



C. The Hopf symmetry description of achiral 
non-Abelian nematics 

1. General aspects 

We haven't shown that all properties of electric, mag- 
netic and dyonic excitations are captured by F(Hm) x 
CHei- For example, we would like F{Hm) x CiJei to 
reproduce the fusion rules of these modes. This can 
be done, by introducing the coproduct A which in turn 
determines the tensor products of the irreps. This 
works as follows: A is a map from F{Hm) x CH^i to 
F{Hrn) X CHel (E" F{Hm) X CH^i, that respects the mul- 
tiplication (i.e. it's an algebra morphism): 

AiFhgPh'g') = AiPhg)AiPH'g') (13) 

Given an element Phg, the coproduct can be written out 
in a basis of F{H„i) x CHei F{Hra) x CHeV- 

^{Ph9)= X! >^}n^giJi2,g2Phi9i<» Ph292- 

/ii,/i2G-ff,„ ,31 ,g2&H^i 

Because this is rather cumbersome notation, we adopt 
the more convenient Sweedler's notation instead. For 
any a S F{Hm) x CiJg; we write 

A(a) = ^a(^) (g)a(2). (14) 

(a) 

This means we can write A (a) as a sum of elements of 
the form a^'^^ ® a'^'^\ with a^^) and a'^) g F(7J„) x CHei- 
Now if we have two irreps 11^ and 11^ , their tensor prod- 
uct is a representation of F{H„i) x CH^i whose action on 
a is given by 

(n^®n^)a^ (n^®nf)A(a) = ^n^(a(i))®nf(a(2)). 

(a) 

It is possible to choose the coproduct in such a way that 
it produces the right fusion rules. The answer is^i^ 

MPhg)^ J2 Phh'~^9®Ph-g- (15) 

F{Hra) X CHel IS also a Hopf algebra- This means there 
is even more structure on F{Hm) x Ci?ei than we have 
defined until now. Here we will only introduce the struc- 
tures that are relevant for this chapter and the next. 
A Hopf algebra has a counit e, which corresponds to 
the trivial or vacuum representation- For the case of 
F{Hra) X CHel, £ is defined by 

e{Phg) = 5h,e- (16) 

One may also consider e as a one-dimensional representa- 
tion of the double F{Hm) x CHei, whose tensor product 
with any irrep 11;^ gives H^: 

n^(g)e~e(g)H^ ~n^. (17) 
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We introduce one more structure: the antipode S', defined 
for F{H,„) X CH,i by 

SiPh9)^Pg-i.h-i9-'. (18) 

It is used to define the conjugate or antiparticle repre- 
sentation of 11;^: 

nj(p,5) = (n^5(F^5))*, (19) 

where t denotes the transpose. The properties of the an- 
tipode imply that this 11^ is a representation, and that 
the vacuum representation e appears in the decomposi- 
tion of nj® n^: 

nj®n^ = £©0nf. (20) 

B,f3 

This property explains the term "antiparticle irrep" : an 
irrep and its anti-irrep can "annihilate" into the vacuum 
representation s (i.e. there is no topological obstruction 
to such a decay) . This discussion applies to general Hopf 
algebras, and therefore to any physical system character- 
ized by a Hopf symmetry. 

2. Braiding and quasitriangular Hopf algebras 

Braiding plays a crucial role in the double symmetry 
breaking formalism. We have discussed braiding as it fea- 
tures in the present context in some detail elsewherei^, 
therefore we will be brief here. We are especially inter- 
ested in the way braiding is implemented in the alge- 
braic structure of a modified quantum double. First we 
review the case where the Hopf symmetry ^ is a quan- 
tum double D(H), which is a modified quantum double 
with Hei — Hjn — H ■ Then we address the case where 
A = F{H„,) X CHei, with Hm + Hei. 

Braiding addresses the following question: What hap- 
pens to a two-particle state when one excitation is adia- 
batically (i.e. slowly) transported around the other? The 
braiding properties are encoded in the braid operator TZ. 
When two defects \g > and \h > {g,h £ Hm — H) are 
braided, it is known what the outcome is (it follows from 
homotopy theory^^-^^'^^). If \g > lives in and \h > in 
V^, then 7^ is a map from (g) to (g) whose 
action is defined by 

11 ■ \g > (»\h >= \ghg~'^ >(^\g > (21) 

Note that it braids the defect to the right halfway around 
the other defect, and we call this a half -braiding. To 
achieve a full braiding, or monodromy we have to apply 

The equation for the braiding of defects \g > and \h > 
we have just discussed applies equally well to the case of 
global as to the case of gauged symmetry. 

Electric modes braid trivially with each other^«: 



The (full) braiding of an electric mode with a topologi- 
cal defect leads to the phase factor causing the famous 
Aharonov-Bohm effect. In the present non-Abelian con- 
text that means that if we carry a particle in a state \v > 
of a representation a of the group H adiabatically around 
a defect with topological charge g E H then that corre- 
sponds to acting with the element g in the representation 
a on \v >: 

n\g > (»\v >^ a{g)\v > (^\g > . (23) 

In the global case, during the parallel transport the par- 
ticle is following a curved path in its internal space. It 
is being "frame dragged", as it is calledSi. To be spe- 
cific, one defines a local coordinate frame somewhere at 
the start of the path in G characterizing the defect. Then 
one lets the elements of the path act on this initial frame, 
to obtain a new frame everywhere on the path. An elec- 
tric mode is then parallel transported around the defect 
if its coordinates are constant with respect to the local 
frames. This is basically the reason that one obtains the 
same outcome for parallel transport as in a local gauge 
theory, the only difference being that in a gauged theory 
the particular path in Hilbcrt space taken by the electric 
mode is gauge dependent and therefore not a physical 
observable. Only its topological winding number leads 
to an observable effect. 

There is a continuum formulation of lattice defects 
in terms of curvature and torsion sources in Riemann- 
Cartan geometrjiS^iSi. In this geometrical approach one 
can also explicitly evaluate the outcome of parallel trans- 
port of an electric mode around a defect. This idea has 
been applied to quite a few phases, such as superfluid 
helium, where the symmetry is also globaP^. It has also 
been applied to uniaxial nematic liquid crystals in the one 
constant approximation's (in the absence of diffusion) . 

One of the advantages of introducing the Hopf symme- 
try A (which we take to be a quantum double for now) 
is that a Hopf algebra is naturally endowed with a so 
called universal R matrix R. R is an element of A(E) A. 
It encodes the braiding of states in irreps of A: to braid 
two states, |0i > in Hi and > in H2, act with R on 
\(j)i > ^\(j)2 >, and then apply the flip operator r This 
composition is called the braid operator TZ: 

7^(|(/)l > «)|02 >) = r o (Hi «) H2) o i? o > (g)\(j>2 > . 

where the action of t is just to flip any two vectors \4>i > 
and \(j>2 >: 

r(|0i >(g)|(/.2 >) = 102 > > • (24) 

If > is in the vector space Vi, and |02 > in V2, then 
\4>i > (8'|02 > is a vector in Vi (g) ¥2- Then 7?.(|0i > 
(g 1 02 > ) is a vector in V2 (g Vi . 

The universal R matrix is an invertible element of ^(g 
A, i.e. there is an R~^ G A^ A which satisfies 



TZ\vi > IS)\V2 >= \V2 > <»\vi > 



(22) 



RR~^ = R^^R = 1 «) 1. 



(25) 
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72, corresponds to braiding the particle on the right in a 
counterclockwise fashion halfway around the particle on 
the left. Using we can define the inverse braiding, 

which is the clockwise braiding of the particle on the right 
halfway around the particle on the left: 

7^-l=i^-loT. (26) 

It is sometimes convenient to write R in Sweedler's no- 
tation: 

E = ^i?(i)0i?(2). (27) 

[R) 

We can let R act on n-particle states. To do this, We 
define 

Rij^^l®---® R^^^ R^^^ • • • 1 (28) 

{R) 

where i?*^^^ is in the i-th, and i?^^) in the j-th position. 
Rij implements the half- braiding of particles i and j. i 
needn't be smaller than j. For example, on a two particle 
state i?2i = E(it) R^^^ ® R''^''- 

For the D{H) = F{H) x CH case, the universal R 
matrix is given by 

R=^Pge(g)g. (29) 

g&G 

The braid operator TZ that is derived from this R repro- 
duces the braiding of the different modes discussed in this 
section. 

The universal R matrix satisfies the Yang-Baxter equa- 
tion: 

R12R13R23 — -^23-^13^12- (30) 




FIG. 1: The Yang-Baxter equation. 



A Hopf algebra with a universal R matrix is called a 
quasitriangular Hopf algebra, and the quantum double 
D{H) is quasitriangular. 

So far we have discussed braiding for the quantum dou- 
ble of a discrete group, but we are also interested in the 
case where Hm ^ H^i, and then we need to know what 
the outcome is of braiding a vector \v > around a defect 
\g > (the braiding of defects is the same as above, see 
(|21|l ). In the case of non-Abelian nematics (and many 



other cases), the vector |t; > is acted on by some element 
oi Hei, and this element is independent of the vector \v >. 
In other words, there is a map F : H„i 1-^ H^i that sends 
g G Hm to F(g), which acts on \v > when |t; > is parallel 
transported around \g >■ F is a group homomorphism, 
and is dictated by the physics of the system we are con- 
sidering. In the cases we have studied (in particular the 
cases relevant for this article), F also satisfies the follow- 
ing relations: 

yg,9i,g2 & H,n,h £ Hei : (31) 
r(5i) ■ 92 = 919291 ^ (32) 
r{h ■ g) = hr{g)h-\ (33) 

F{Hm) X CHel is then a quasitriangular Hopf algebra 
with the following braid matrix: 

R= ^9e®F(g). (34) 

This equation precisely encodes what we have described 
above. The inverse of R is 

R-' = Y,Pge(g>Tig-'). (35) 

The quantum double D{H) is a special case of a 
generalized quantum double, with H,n = H^i = H, 
h ■ g = hgh~^ yh,g g H, and F = id, the identity 
operator. We will see examples of phases with nontrivial 
F later on. 



3. Achiral non-Abelian nematics 

We will now explicitly describe the Hopf symmetry 
relevant for non-Abelian nematics with tetrahedral, 
octahedral and icosahedral residual symmetry. The 
Hopf symmetry A is of the form discussed above: 
A — F{Hm) X CHei- We will exphcitly describe the 
electric and magnetic groups, and the defect classes. 
We also analyze the consequence of the presence of 
inversions and reflections in H^i- 

Achiral tetrahedral nematic 

The electric group is H^i — Td, the group of sym- 
metries of a tetrahedron, including reflections (since the 
phase is achiral), see figure [3 We denote elements of Td 
as permutations of the four vertices of a tetrahedron, e.g. 
(12), (134), (13)(24), etc. 

Before we can define the magnetic group, we first de- 
scribe a common parameterization of SU{2), and the 
two-to-one homomorphism from SU{2) to 5*0(3). To 
specify a rotation in SO(3), one specifies an axis around 
which the rotation takes place, and a rotation angle 
— TT < 6* < vr). Denote by n a unit vector along the 
axis of rotation, and define positive 6 to correspond to 
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FIG. 2: The tetrahedron, with labelled vertices. 



counterclockwise rotation with respect to fi. Then this 
rotation is denoted by R(fi,9). In this parameterization 
we can envisage SO{i) as a ball of radius tt with antipo- 
dal points on the surface identified. 

We can parameterize matrices in SU (2) in a very 
similar way: take any unit vector n, and any angle 
— 27r < 6 < 2tt). Notice how 9 runs over a larger range 
than in the 50(3) case, we now have a sphere with radius 
2tt and the surface of the sphere corresponds to the cen- 
ter element —1). Now associate to h and 9 the following 
matrix in SU (2): 



Preferred element qa 


Conjugacy class A 


ite 


(^±e — ±te| 


±[123] 


C±[123] ~ 

-i-/rioQi ri 9/11 ri/ioi ro/iQi\ 


±[124] 


C±[124] ~ 

±{[124], [132], [234], [143]} 


[(12)(34)] 


C[(12){34)] = 

{±[(12)(34)],±[(13)(24)], 
±[(14)(23)]} 



TABLE I: The conjugacy classes of T, and the preferred ele- 
ments gA- When we write ± in front of a class, we mean that 
there are two classes, one with plus signs in front of all the el- 
ements, and one with minus signs in front of all the elements. 
On the other hand, if we write ± in front of an element inside 
a class, we mean that the both the element with a plus and 
minus sign in front of it are present in the class. 



u{n^ 9) — exp{i—h ■ a) — cos{ — )l + isin{—)n ■ a (36) 

where a = (ux, fj,, cr^) is a vector containing the three 
Pauli matrices, and 1 is the unit matrix. 



1 

1 l^'^ y- 



-i 

1 



1 

-1 



(37) 

5*0(3) is now easily 



The homomorphism p : SU (2) 
accomplished: replace u by R. 

In the achiral tetrahedral phase, the magnetic group is 

f: 

ni(0(3)/Td) = ni(50(3)/r) - ni(5c/(2)/f ) ~ f. 

T is the inverse under p of T, the tetrahedral group in 
50(3). Global symmetry transformations act on the de- 
fects in T, so the defects are grouped together in orbits 
under the action of the electric group, in this case T^. We 
note that this is not the same as the conjugacy classes 
of T\ To see this, we must first fix a notation for the 
defects. Write an element of T as 



u{h,9) = sgn(cos—)u' (n, 9') 



-TT <9' <TT 

9' = 9{mod2Ti). (38) 



Every defect corresponds to a u(n, 9) (the endpoint of 
the path in SU{2) that characterizes the defect). We 
denote the defects as cycles of T with square brackets 
with a plus or minus sign, e.g. ±[123]. The minus sign 
corresponds to the 2tt defect, i.e. the nontrivial loop in 



50(3). The defect [123] corresponds to the u'(n, 9') with 
—TT<9' < TT such that 

R{n,9') = (123). 

We need to know the axis ni in figure|21 such that (123) 
corresponds to a ^ rotation around rii. 

From the figure we see that lii = -^(1,1,1). Thus we 

define [123] — u{ni,^). The trivial defect is denoted by 
1, and the 2tt defect by —1. 

To have a notation for all the defects in T, we first 
define the following axes: 



ni 



n-3 



n2 



ri4 



Then the defects are given by ±e and 



± [123] 


— ±u{ni, 


27r 


27r 

±[124]=±^n2,-) 


±[124] 


= ±u{ri2, 


27r 


2tt 

±[234] -±w(n3,y) 


±[134] 


— ±u(n4, 


27r 


±[132] =±M(ni,-y) 


±[142] = 


±u(?f2, - 


2tt 
'J' 


±[243] =±M(n3,-y) 


±[143] = 


±u(?f4, — 


2tt 


±[(12)(34)]=±u(z,^) 


±[(13)(24)]=±u(y,7r) 


±[(14)(23)]=±u(i;,7r) 
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Classes of T under the action of Td 



C-e 
C[1231 ^ C'[124] 
C-[1231 ^ ^"-[124] 
C[(12)(34)1 



TABLE II: The defect classes in a theory with Td symme- 
try. The classes are unions of conjugacy classes of T. The 
conjugacy classes of T are defined in table |T| 



Using our notation, we can determine the defect classes 
under the action of Td- First consider an element 
R{n,(l)) of T. Its action on a defect u{ni,9) gives 
u{n, (j))u{ni, 9)u{n, 4>)~^- 

Now consider transformations in Td that are not con- 
nected to the identity, such as the element (12). These 
are called "large" symmetry transformations. We can 
always write a large symmetry transformation as Inv x 
R{n, 0), where Inv is the inversion matrix: 



Inv 











(39) 



Inv acts trivially on all the defects, thus the ac- 
tion of Inv X R(h,(j)) on a defect u{ni,9) gives 
u{n, <p)u{ni, 9)u{n, <p)~^. 

We will now write the large symmetry transformations 
as Inv X R{h, (j)). First we define the following axes: 



mi 



TO5 



1 

1 

"71 
1 



(1,-1,0) m2 = 
(0,1,1) m4 = 

(1,0,1) TO6 



1 

71 
1 

71 
1 



(-1,0,1) 
(0,-1,1) 
(1,1,0). 



' ° 
The inversions arc given by 

(12) = Inv X R{rhi,n) (13) = Inv x R{m2, tt) 
(14) = Inv X i?(TO3, tt) (23) = Inv x i?(m4, tt) 
(24) = Inv X R{rh5, tt) (34) = Inv x R{ihQ, tt). 

We can now derive the multiplication table of T using 
the multiplication in SU{2). For example, 

[123] [123] = -[132] 
[123] [124] = [(14)(23)] 
(12)[(12)(34)](12) = -[(12)(34)] 
etc. 



yo 


C6iitr£Lliz6r Ncy 


IrrGps 


±e 


T 




±[123] 


{e,(123),(132)} -Za 




±[124] 


{e,(124),(142)}~Z3 




[(12)(34)] 


{e,(12)(34),(13)(24),(14)(23)} 

~ Z4 


n±[^^^fceZ4 



TABLE III: The irreps of F{f ) x CTd. 
Td. 



The a are irreps of 



The Hopf symmetry of the achiral tetrahedral nematic 
is the modified quantum double F{T) x CTd- We just de- 
fined the action of Td on T, which sets the multiplicative 
structure of F(f') x CTd. 

Now that we have elucidated the action of Hei on Hm, 
we can determine the defect classes. These defect classes 
are the union of conjugacy classes of T. The conju- 
gacy classes of T are shown in table and the defect 
classes in table m The "small" symmetry transforma- 
tions (that are connected to the identity) simply conju- 
gate the defects, while the large symmetry transforma- 
tions may transform defects in different conjugacy classes 
into each other. 

The centralizcrs and irreps of F{T) x CTd are given in 
table Iml 

We must carefully interpret the result that defect 
classes can be larger than conjugacy classes. The de- 
fects in the same class have the same energy, since the 
symmetry transformations commute with the Hamilto- 
nian. They also have isomorphic centralizcrs. However, 
the cores of defects only related by a large symmetry can- 
not be interchanged. This is due to the fact that the 
symmetry isn't connected to the identity. Thus these 
defects are not topologically equivalent: they cannot be 
deformed into each other with a finite amount of energy. 
So we might conclude that we shouldn't act with large 
symmetries on the defects, and we should work with con- 
jugacy classes. However, we cannot simply neglect the 
large global symmetries, since they act on the electric 
modes. 

Finally, we need a map F from Hm to Hei to define a 
braiding (see above): F is given by 

F(±[123]) = (123),F(±[(12)(34)]) = (12)(34), etc. 

Thus F turns square brackets into round brackets and 
neglects the eventual minus sign. 
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Pref. el. gA 


Defect class A of O 




0±e — \±e| 


±[123] 


C±[i23l = ±{[123], [124], [142], [132], [134], 
[234J, [i42J, [i43J, [243J| 


It [lZo4J 


[1342], [1423], [1432]} 


[(12)(34)] 


q(i2)(34)i = {±[(12)(34)],±[(13)(24)],±[(14)(23)]} 


[12] 


C[i2i = {±[12], ±[13], ±[23], ±[14], ±[24], ±[34]} 



TABLE IV: The defect classes of O under the action Oi, and 
the preferred elements gA- See |l] for a description of the no- 
tation. 



Achiral octahedral nematic 

The electric group is Oi, which is the octahedral group 
O (consisting of the rotational symmetries of a cube), 
plus the inversions and reflections of a cube. O is iso- 
morphic to 5*4 (all permutations of the diagonals of a 
cube). Thus we can write elements of O as cycles, e.g. 
(1234), (123), (23), (13)(24), etc. Defining Inv as in equa- 
tion (|39|l . we have that Oi — {O, InvxO}. Thus elements 
of are (1234), (123), Inv, Inv x (134), etc. 

In analogy with the achiral tetrahedral nematic dis- 
cussed above, we can denote elements of the magnetic 
group O as cycles with square brackets, with an eventual 
minus sign. Examples are ±[123], ±[(12)(34)], ±[(12)], 
where the minus sign is the 27r defect. 

The Hopf symmetry is F{0) x COi. The defect classes 
in O under the action of Oi are given in table IIVI The 
inversion Inv acts trivially on the defects, so for the octa- 
hedral nematic the defect classes are in fact the conjugacy 
classes of O. In general, when a group carries a sublabel 
i, it means that the group contains Inv. 



Pr. el. QA 


Defect class A. of I 


±e 


C±e = ±{e} 


It [l^oj 


/"* . , -i-/rio'^i ri "^o! ri 0/11 ri/i9i riof^i 

'-'±[1231 ^ ±ll-'-^'jJi [J-'J^Ji [-1-^4], [i42J, [i^Oj, 

[152], [134], [143], [135], [153], 
[145], [154], [234], [243], [235], 
[253], [245], [254], [345], [354]} 


±[12345] 


C±[i2345] = ±{[12345], [12453], [12534], [13254], 
[13425], [13542], [14235], [14352], 
[14523], [15243], [15324], [15432]} 


±[12354] 


C±[i2354i = ±{[12354], [12435], [12543], [13245], 
[13452], [13524], [14253], [14325], 
[14532], [15234], [15342], [15423]} 


[(12)(34)] 


C[(12)(34)l = 

±[(13)(25)], ±[(13)(45)], ±[(23)(45)], ±[(14)(23)], 
±[(14)(25)], ±[(14)(35)], ±[(24)(35)], ±[(15)(23)], 
±[(15)(24)],±[(15)(34)],±[(25)(34)]} 


±[123] 


C^±[i23l = ±{[123], [124], [125], [132], [134], 
[234], [135], [235], [142], [143], 
[243], [145], [245], [345], [152], 
[153], [253], [154], [254], [354]} 



TABLE V: The defect classes of / under the action and 
the preferred elements gA- See table|I]for a description of the 
notation. 



the minus sign is the 27r defect. The Hopf symmetry is 
F{I) X Cli- The defect classes in / under the action of li 
are given in tabled The inversion Inv acts trivially on 
the defects, so the defect classes are the conjugacy classes 
of /. 



Achiral icosahedral nematic 

The electric group is li, which consists of the icosa- 
hedral group / (the rotational symmetries of an icosahe- 
dron), plus inversions and reflections of the icosahedron. 
/ is isomorphic to (the even permutations of the five 
cubes inscribed inside an icosahedron). Thus we write 
elements of / as cycles, e.g. (12345), (14)(23), etc. Just 
as in the octahedral case, — {I, Inv x /}. 

We can denote elements of the magnetic group / as 
cycles with square brackets, with an eventual minus 
sign. Examples are ±[123], ±[(12) (345)], ±[(15)], where 



III. DEFECT CONDENSATES AND HOPF 
SYMMETRY BREAKING 

A. The Hopf symmetry breaking formalism 

Using the Hopf symmetry description of a phase, one 
can study phase transitions induced by the condensation 
of any mode, be it electric, magnetic or dyonic. The 
theory of Hopf symmetry breaking was first proposed 
by Bais, Schroers and Slingerlandii^. who developed 
the framework that determines the physics of the broken 
phase, and applied it to discrete gauge theories. The aim 



11 



of the present article is to apply this framework to non- 
Abelian nematics. It turns out that we were physically 
motivated to alter one step of the framework, namely the 
definition of the residual symmetry algebra. In this sec- 
tion, we briefly introduce the Hopf symmetry breaking 
formalism, relying strongly on physical motivation. The 
mathematics of our approach, which differ slightly from 
the mathematics in papers just mentioned, are discussed 
in another paper—. 

Classically, a condensate corresponds to a nonzero 
expectation of some (dis)order parameter field. This 
nonzero expectation value has certain symmetries, its 
symmetry group is called the residual symmetry group 
H. The quantum interpretation is that the ground state 
of the system corresponds to a non-vanishing homoge- 
neous density of particles in some particular state. The 
residual symmetry operators are the operators that leave 
that particular state invariant. Whether we are working 
classically or quantum mechanically, the residual symme- 
try operators are determined the same way: determine 
the operators that act trivially on the condensate, which 
is a vector in some irrep of the original symmetry. 

If the original symmetry was a group G, it would be 
easy to define what we mean by the symmetry operators 
of a vector |(/)o > in an irrep 11 of G, namely operators 
h E G that satisfy 



h- 100 \(t)o > 



(40) 



If the original symmetry is a Hopf algebra A, on the other 
hand, the definition of "residual symmetry operators" is 
not so trivial. This is discussed in a related paper'^: There 
is a physically very attractive definition of a "residual 
symmetry operator" : an operator is a residual symmetry 
operator if its action on a particle state is not affected by 
the fusion of that particle with the state of the particles in 
the condensate. This means that if this operator acts on 
any particle, and that particle fuses with the condensate, 
then the action of the operator on the particle before 
and after fusion with the condensate is the same. We 
have to make a choice at this point: whether the particle 
fuses with the condensate from the left or the right. This 
is an important point as we argue elsewherei*. Here we 
choose fusion with the condensate from the right. This 
leads to the definition of the right residual symmetry 
algebra, which consists of all symmetry operators in A 
whose actions on a state are not affected by fusion of the 
state with the condensate on the right. %■ is the analog 
of the residual symmetry group when a group symmetry 
is spontaneously broken. We note that there are cases 
where %■ isn't a Hopf algebra. 

Once %. has been established, we must take a close look 
at the particles in this broken theory, which we now con- 
sider to be the irreps of 7^ . It turns out that some parti- 
cle species don't braid trivially with the condensate. The 
presence of such particles in the system implies a half-line 
discontinuity in the condensate, across which the internal 
state of the condensate jumps, which means that these 
particles have to be connected to a domain wall. Hence, 



Single defect condensate in 
A = -F'(f') X CTd 


K 


% 


U 


\-e> 


Ci 


F(T) X CTd 


F(T) X CTd 


![123] > 


Cs 


F{T/C3) X CCa 


Die) 


|[(12)(34)] > 


C2 


F(T/C2) X CD2 


D[C2) 


1 - [(12)(34)] > 


C2 


F(T/C2) X CDa 





TABLE VI: Single defect condensates in a tetrahedral ne- 
matic. A is the original Hopf symmetry, K is the smallest 
subgroup of Hm that contains all the defects gi that are in 
the condensate, % is the residual symmetry algebra, and U 
is the unconfined Hopf symmetry. 



such particles will be confined. We can similarly deter- 
mine which particles are unconfined, i.e. braid trivially 
with the condensate, and (in the cases we've worked out) 
these particles turn out to be the irreps of a new Hopf al- 
gebra called the unconfined symmetry algebra which will 
be denoted as U. 

Thus, in contrast with the conventional symmetry 
breaking analysis, we have to distinguish two steps in 
the symmetry breaking: first to the residual symmetry 
algebra 7^, and then to the unconfined symmetry algebra 
U. One should not make the mistake of believing that 
U contains all the information about the broken phase: 
one might believe that the confined particles should be 
neglected because they cost an infinite amount of en- 
ergy to create in a system of infinite extent, since the 
half-line discontinuity costs a finite amount of energy per 
unit length. However, the half-line may end on another 
confined excitation, thus giving rise to a wall of finite 
length. We call a configuration consisting of confined 
excitations connected by walls, such that the overall con- 
figuration is unconfined, a hadronic composite, in anal- 
ogy with hadrons in Quantum Chromodynamics, where 
hadrons are unconfined composites of confined quarks. %■ 
contains all the information about the hadrons, although 
it is not trivial to extract this informatioFj^ experimen- 
tally 



B. Defect condensates and residual symmetry 
algebras 

One can derive general formulae for %■ and W in the 
case of defect condensates in a phase with F{Hm) x CiJg; 
symmetry^. Here we will just summarize the resulting 
formulae, and in the next section apply them to the classi- 
fication of all possible defect condensates in achiral tetra- 
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Class sum def. condensate 
of ^ = F{T) X CTd 


K 


Tr 


U 






r yi ) )< d 


r (1 ) X lUJ d 


\C[12S\ >, \C-[123] >, 


f 


CT 


D{e) 


\C[124] >, |C'_[i241 > 








|C[(12){34)] > 


D2 


F{Zs) X CTd 


F{Z:i) X C(Z3)d 



TABLE VII: Conjugacy class sum defect condensates in an 
achiral tetrahedral nematic. (Zs)^ is isomorphic to the per- 
mutation group of 3 elements. 



hedral, octahedral and icosahedral nematics. 

Given a phase with Hopf symmetry A = F{Hm) x 
CHei, the condensate |<^o > is a vector in some irrep 11;^ 
of A. We must demand that {(jjo > braid trivially with 
itself, otherwise the condensate is ill-defined (the con- 
densate itself would be filled with half- line singularities). 
Thus we demand that 

7^ • 100 > |</.o >= \<f>o > I'/'o > • (41) 

Now let us take |0o > to be a defect condensate, which 
means it is a vector in a magnetic irrep (1 is the 
trivial irrep of the normalizer Na)- There are different 
types of defect condensates which we wish to analyze. A 
basis of the vector space on which this irrep acts is given 
by {\gf >}, where the gf are the different defects in A. 
We consider the following types of condensates: 

• Single defect condensate 

Igf > (42) 

• Class sum defect condensate 

100 >- J2 \9^ >=■■ > (43) 

where Cg^ is a defect class (i.e. the orbit of gA 
under the action of Hei). We denote the condensate 
by |Cg^ >, where gA is the preferred element of A. 

• Combined defect condensate 

I'/'o >= E \9^ > (44) 

where £' is a subset of the defects in one class. We 
need only take the elements to be within one class 
because it turns out that we need only study the 
cases where the condensate is the sum of vectors in 
the same irrep (all other phase transitions can be 
interpreted as a sequential condensation of vectors 
in different irrepsi^.). 



Combined def. cond. 
of ^ = F{f) X CTd 


K 


% 


U 


1[123] > +|[132] > 


Ci 


F{T/C3) X CC3„ 


CCit, 


|[(12)(34)] >+|[(13)(24)] > 


D2 


F(T/D2) X CC2v 


CCi. 


{|[(12)(34)] > + 
|[-(12)(34)] >+|[(13)(24)] > 

1 1 Q^^'0/!M \1 

+ 1 - l(.lj)(^4)J >| 


D2 


F{T/D2) X CD2d 


CCiu 


|[(12)(34)] >+|h(12)(34)] > 


C2 


F{T/C2) X CD2d 


F{C2) 

XCC2v 


|[(12)(34)] >+|[(13)(24)] > 
+ |[(14)(23)] > 


D2 


F{1z) X CT 


D{13) 



TABLE VIII: Combined defect condensates in an achiral 
tetrahedral nematic (that satisfy trivial self braiding). No 
two defect condensates have simultaneously the same Tr and 
U. Thus in principle the different defect condensates are dis- 
tinguishable. 



The single defect and class sum defect condensates are 
a special case of combined defect condensate. 



1. Single defect condensate 

Condense \gA > in the magnetic irrep 11^. The con- 
densate \gA > satisfies the trivial self braiding condition 
mentioned above. The residual symmetry algebra is 

Tr = F{H„J{gA))®CNA (45) 

where we define (gA) to be the smallest subgroup of Hm 
that contains gA- 

This result for 7^ has a very natural interpretation: 
the residual electric group is Na, the subgroup of H that 
doesn't conjugate the defect. The magnetic part H/{gA) 
is not necessarily a group. It consists of left cosets of 
(qa) = {• • ■ 9a^: e, gA, g\, ■ ■ •}■ The defects are now de- 
fined modulo the condensate defect \gA >■ In other 
words, if a particle in a magnetic irrep of the residual 
symmetry % fuses with the condensate \gA >, it is left 
unchanged. Thus its topological charge is defined modulo 
9A- 

One can prove that 7^ is a Hopf algebra <;=^ (gA) is 
a normal subgroup of Hm Hm/igA) is a group. 

The unconfined symmetry algebra is 

U = FiT-\NA)/{gA)) X CiNA/TigA)). (46) 
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Single def. cond. 
oi A = F{d) X CO, 


K 


% 


U 


|-e> 


Ci 


F{0) X CO, 


F{0) X CO, 


1 ± [12] > 


C2 


F{0/C2) X CD2^ 


F{C2) X C02, 


1 ± [123] > 




F{0/C'3) X COs, 


CO, 


± [1234] > 




F{0/C4.) X CCu 


CO, 


|±[(12)(34)] > 


C2 


F{0/C2) X CDi, 


F{D2) X CD2, 



TABLE IX: Single defect condensates in an achiral octahedral 
nematic. When we write ± in front of a condensate, we mean 
that the condensate with or without the minus sign gives the 
same symmetry breaking analysis. 



Note that we condensed \gA >, where gA was a chosen 
defect in the defect class A. gA was chosen arbitrarily, so 
our formulae are general. 



2. Class sum defect condensates 
Condense the sum of the defects in the defect class A: 



\<t>o >= E \9f \Cg.. > ■ 

A class sum defect condensate satisfies the trivial self 
braiding condition (|41|l : 

= 7^( ^ \gf > ® ^ \g^ >) 

gfeA g^eA 

= E(E \9f9ti9fr'>)®\9f> 
gfeA gfeA 

= E ( E \9t >) X \gf > . 

gfeA gfeA 
= \Cgj, > ®\Cg^ > 

In going from the second to the third line, we use the fact 
that gAg^^ — A for any g e Hm- 

A class sum condensate doesn't break the electric 
group at all! Namely, all elements of H^i act trivially 



Class sum def. cond. 
of A = F(d) X COi 


K 




u 


|C[(12)(34)1 > 


D2 


F{Da) X CO, 


F(D3) X CDs, 


0±[i23] > 


f 


^^(02) X CO, 


F(02) X CO2, 


[^±[1234] > 


O 


CO, 


CO, 


lQl21 > 


O 


CO, 


CO, 



TABLE X: Conjugacy class sum defect condensates in an achi- 
ral octahedral nematic. When we write ± is front of a con- 
densate, we mean that the condensate with or without the 
minus sign gives the same symmetry breaking analysis. 



on a defect class, since for any g e Hgi we have 

<? • 100 >= 5 • ( E \9f >) = 

gfeA 

= E If • 5.^ > = E \9f >= > ■ 

gfeA gfeA 

Thus this condensate is invariant under all of Hei- For 
this reason, such a condensate is admissible in a theory 
where the symmetry is gauged and we call it a gauge 
invariant magnetic condensate (the condensate respects 
gauge invariance). In a global theory all condensates are 
admissible. 

The residual and unconfined symmetry algebras are 
respectively 

% = F{H„JK) X CHei (47) 
U = F{H„JK) X CH,i/T{K), (48) 

where K is the smallest subgroup of Hm that contains 
the class A. From this definition, one can prove that K 
is a normal subgroup of Hm- Thus H„i/K is a group, 
and 7^ is a Hopf algebra. 

Later on, we will consider conjugacy class sum deject 
condensates, i.e. condensates of a sum of defects in the 
same conjugacy class of Hm- This can be a whole defect 
class, or it can be smaller than a whole defect class. If it 
is smaller, then the electric group is partially broken. 



3. Combined defect condensates 

Start with a phase with F{H„i) x Ci/e; symmetry. 
Choose an irrep 11^, and consider a condensate of the 
form X)g e_B \9i with E a subset of the defects in one 
defect class. 
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Combined def. cond. 
oi F(d) X COi 


K 


% 


U 


|[123] > +|[132] > 




F{0/C-a) X CAm 




1 n 9i \ -1-1 f'^/ii \ 




r [{J / U2 ) X ^UAi 




|[(12)(34)] > 

+ i[(13)(24)] > 


D2 


F{0/D2) X CDi, 


F(C2) X CC2> 


[12] > +|[13] > +|[23] > 




F{0/Da) X CD3» 





TABLE XI: Combined defect condensates in an achiral octa- 
hedral nematic. 



The demand of trivial self braiding (|41|l gives 
'^{Y.\9^> \9k >) = I] ^ \9k > 

gz^E gkGE g,eE gk&E 

Yl \9^9k9i^ > ®Wi > =YY\9'> > 

gk£Egi£E g,eE g^eE 

y9^'^E : {g,gk97^}g,eE = {9k}g,eE- (49) 

It is interesting in itself to study how many different de- 
fect condensates satisfy this criterion. Defect-antidefect 
condensates \g > +\g~'^ > always satisfy this criterion2&, 
as do any set of commuting elements in a certain conju- 
gacy class, and class sum defect condensates. The trivial 
self braiding condition plays a crucial role in determining 

The results for %. and U are: 

% = F{Hm/K) X CMe 

U = F{Ne/K) X CMe/T{K), 

where we must still define all the notation in these results. 

Define the following subset of H^i (which needn't be a 
subgroup) : 

Ve C Ha : Ve = {x,NA}g^^E (50) 

where Na C i?e; is the normalizer of the chosen preferred 
element gA in A, and the Xi satisfy XigAX~^ = gi € E. 
Ve corresponds to the set of left Na cosets that corre- 
sponds to the defects in the condensate (under a corre- 
spondence discussed above). 

Define the following subgroup of Hei : 

Me C Hei : Me = {m e Hei : {m • gt}g,eE = {ffilsiS-E} 
= {me Hei : mVE = Ve}. (51) 

Me is composed of the global symmetry transformations 
that leave the condensate invariant. 



A = F(/) X Cli 


K 


% 


U 


\-e> 


Ci 


F{I) X CI, 


F(I) X CU 


|[123] > 


C3 


Fil/Cs) X 




1 - [123] > 


C3 


Fii/C-i) X CC3. 




1 ± [(12)(34)] > 


C2 


F(//C2) X CD2^ 


F{C2) X CC2^ 


[12345] > 


C5 


F{I/Cs) X CC5, 


ca 


- [12345] > 




F{i/Cs) X 


ca 



TABLE XII: Single defect condensates in an achiral icosahe- 
dral nematic. 



Also define 

Ne C H„i ■■ Ne = {ne H,n ■■ {ngiU^^jg^i^E = {gtjg.eE}- 

Using 1(23): T{gi) ■ g2 = gi523r^ Vgi,52 e H^, we can 
prove that 

T-\Me)=Ne and T{Ne) ^ Me- (52) 

Finally, we need one more definition: 

K C : K - ({.gjg.es), (53) 

where ({ffi}gig_E) is the smallest subgroup of Hm that 
contains all the gi E E, i.e. the defects in the condensate. 

The trivial self braiding equation (|49|l implies that K C 
Ne- Thus, according to l(S2Jl and 

TiK) C Me. (54) 

Summarizing, the unconfincd magnetic group is 
Ne/K, and the unconfincd electric irreps are those that 
have r{K) in their kernel, which means that the electric 
group is Me/T{K). If we take a quantum double D{H) 
[H = iJej = H„i), the unconfincd symmetry algebra be- 
comes U = D{Ne/ K), because in that case Me = Ne. 

C. Non-Abelian condensates in liquid crystals 

We are now in a position to apply the results obtained 
in the previous sections to the case of non-Abelian nemat- 
ics. We have worked out pretty much exhaustive listings 
of all possible phases characterized by defect condensates 
in achiral tetrahedral, octahedral and icosahedral nemat- 
ics. In the corresponding tables of defect condensates we 
give, A is the original Hopf symmetry, K is the small- 
est subgroup of Hfn that contains all the defects gi that 
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A = F{i) X Cli 


K 




u 


[(^±[12345] > 


i 


Ch 


ca 


|C± [12354] > 


i 


Ch 


ca 


|C[(12){34)] > 


I 


Ch 


ca 



TABLE XIII: Conjugacy class sum defect condensates in an 
achiral icosahedral nematic. 

are in the condensate, % is the residual symmetry alge- 
bra, and U is the unconfined Hopf symmetry. The defect 
condensates satisfy trivial self braiding, as we required 
above. 

By looking at all the tables of defect condensates in 
this section, we note that two different condensates in 
our tables never give simultaneously the same 7^ and U. 
Some condensates give the same unconfined symmetry 
algebra, but %. is then different^i. Thus there are dif- 
ferences in the spectrum, though these are often hidden 
in the unconfined spectrum of the hadrons corresponding 
to different condensates. These differences are therefore 
quite subtle and may be hard to detect but in principle 
they are distinguishable. The problem with measuring 
defect condensates, for example, is that the conventional 
measuring techniques can measure the electric symmetry 
group (by looking at Bragg reflections, for example), but 
as far as we know there are no techniques yet to measure 
the magnetic symmetry group. Naively this amounts to 
identifying the surviving unconfined degrees of freedom 
and their interactions, for example by certain interference 
experiments. We would need to measure non-Abelian 
statistics to probe the braiding properties of the parti- 
cles in the broken phase. Only recently have there been 
direct measurements of fractional statistics21. Neverthe- 
less, if suitable techniques were developed, then we could 
use our tables to identify a plethora of new phases and 
determine which condensates they correspond to. 



A = F{i) X Cli 


K 


% 


u 






-T (_i 1 JJ2 ) X lU-L/'2i 




|[(12)(34)] + |[(13)(24)] > + 
_i_i r/i /I ^ 




F{I/D2) X CT, 


CCs, 


|[(12)(34)] >+|[(12)(35)] > + 
+ |[(12)(45)] > 




Fil/Ds) X CDs, 


ca 


{|[(12)(34)] > +|[(13)(25)] > + 
-|-|[(i5)(24)J > -|-|[(23j(45)J > 
+ |[(14)(35)] >} 
|[123] > +\[132] > 


Db 

L^3 


F{I/Ds) X CDs, 

r /O3 ) X 'i^Uzi 


ca2. 


|[123] > -h|[134] > + 
+ |[142] > +\[24S\ > 


f 


F{I/T) X CT, 


ca 


+ |[132] > +|[134] > + 
+ |[234] > +|[142] > + 
+ |[143] > +|[243] >} 


T' 
1 


TT'f T l'T'\ (T^T' 




||[lzoJ > -t-|[10zj > -h 

+ |[135] > +|[253] > + 
+ |[142] > +|[134] > + 
+ |[243] >} 


f 
1 






{| [12345] > -)-! [13524] > 
+ [[14253] > +1 [15432] >} 




F(I/C5) X CDs^ 


ca. 



TABLE XIV: Combined defect condensates in an achiral 
icosahedral nematic. 



1. Achiral tetrahedral nematic 

We have listed all defect condensates in an achiral 
tetrahedral nematic: the single defect condensates are 
collected in table IVII the class sum condensates in table 
IVIII and the combined defect condensates in table I VI III . 



2. Achiral octahedral nematic 

The single defect condensates breaking F{0) x COi 
are given in table Hxl The class sum defect condensates 



are given in table |X] Finally, the combined defect con- 
densates are given in table IXII We note that the list 
presented here is very representative. The other conceiv- 
able defect condensates give the same %. and U as one 
of the defect condensates shown here (except for a small 
difference: there may be condensates where K is actually 
the double of a given here. That slightly changes the 
magnetic part of %, but doesn't affect Li). These other 
defect condensates are trivially different from the ones in 
the table: for example, they may be permutations of the 
of the numbers used in the naming of the defects. 
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3. Achiral icosahedral nematic 

The single defect condensates in F{I) x C/,; are given 
in table IXIll The class sum defect condensates are given 
in table Ixllll Finally, the combined defect condensates 
are given in table IXIVI The sample we present in the 
table is basically exhaustive. 

4. Comments on the conjectured phases 

The tables given above, containing the possible phases 
induced by defect condensation yield a lot of informa- 
tion of both physical and mathematical interest. Yet, 
they don't tell the full story, as they do not describe the 
hadronic composites which we alluded to earlier. Infor- 
mation on these composites is laborious but straightfor- 
ward to extract, and we would like to comment on what 
kind of analysis would reveal the hadronic content of a 
particular phase. 

First we note that more often than not, the residual 
symmetry algebra %. is not a Hopf algebraiS. This is 
perfectly acceptable, since there is no physical reason to 
assume that 7^ is a Hopf algebra. The reason is subtle^: 
for a symmetry algebra to be a Hopf algebra, it is neces- 
sary that the fusion of particles (i.e. the tensor product 
of irreducible representations) be associative. Now some 
of the irreducible representations of %■ correspond to con- 
fined excitations, which means that the condensate is in a 
different internal state to the left and right of the excita- 
tion. Thus particles to the left of this confined excitation 
"see" a different condensate. This leads to the necessity 
of introducing an ordering when taking the tensor prod- 
uct of representations, which corresponds to specifying 
in which order the particles are brought into the system. 

The unconfined symmetry algebra U corresponds to 
the symmetry algebra whose irreducible representations 
are precisely the unconfined representations of 7^. These 
unconfined representations don't suffer from the necessity 
of introducing an order, and therefore we expect U to be 
a Hopf algebra. In all the cases we've worked out this 
turns out to be the case. 

From the tables we learn that different defect con- 
densates may induce the same U, however, they do 
lead to a different {Tr,U) pair. If two phases have 



the same their low energy degrees of freedom share 
many properties (e.g. their representation theory, their 
braiding properties). So to tell these two phases apart, 
it may be necessary to probe unconfined composites 
of confined excitations (which we called hadronic com- 
posites). These may occur at a higher energy scale 
(depending on the precise dynamics). The constituent 
structure of the hadronic composites can be derived from 
%. To determine the admissible composites, one must 
take tensor products of several confined excitations and 
decompose the product into a direct sum of irreducible 
representations. Every unconfined representation that 
appears in such a decomposition, corresponds to a 
hadronic composite. Note that even ii U = Ce, the 
trivial Hopf algebra, there can still be nontrivial com- 
posites, the decomposition then has to yield the trivial 
representation (of %). It is rather tedious to work 
out the sets of simplest allowed composites, because 
the calculation is complicated by the fact that 7^ is 
not a Hopf algebra, but it is straightforward. We have 
refrained from carrying out such an analysis at this stage. 



IV. CONCLUSIONS 

In this paper we have applied the formalism for break- 
ing quantum double symmetries by defect condensates to 
some classes of rather exotic non-Abelian nematics. We 
found a wide variety of conceivable phases each charac- 
terized by a set of unconfined degrees of freedom associ- 
ated with an unconfined algebra U. There may also be 
confined degrees of freedom described by an intermediate 
symmetry algebra 7^. Clearly, whether such phases will 
actually be realized in nature depends on the detailed 
dynamics of these media. It would of course be of great 
value to look for experimental parameters by which these 
phases could be induced and furthermore to develop ob- 
servable signatures by which they could be distinguished. 
These important questions deserve serious attention but 
are beyond the kinematical scope of this paper. 
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^® We give a tailor made summary of these double algebras 
in Appendix A of a related papes-. 

As our treatment is quantum mechanical, the scalars of the 
vector space are in C. 

In the mathematical literature this product is often writ- 
ten as a "bow tie", F{Hm) ixi CHei indicating that there is 
a nontrivial action defined of the factors on each other. 
Formally the structure corresponds to an example of a 
bicrossproduc*— of the Hopf algebras F{Hm) and CHei. 
To keep the notation simple, we use an ordinary product 
sign, to clearly distinguish it from the tensor product sign 
which we will use to describe multiparticle states. 
We actually need to define Td, because Td is a subgroup 
of 0(3), not of 50(3). We get double covers of subgroups 
of SO{3) using the canonical two-to-one homomorphism 
from SU{2) to SO(3). Exploiting the fact that 0(3) ~ 
SO{3) X Z2, we can lift Td to a subgroup of 517(2) x Z2. 
The normalizer of qa is the subgroup of Hei whose elements 
h satisfy h ■ qa ~ QA. 

This follows because F{Hm) x CHei is a transformation 
group algebra. We can then use a theorem described 
elsewhere-. 

Actually, if we are braiding two indistinguishable electric 
particles, then the wave function of the system picks up a 
phase factor e^'^^", where s is the spin of the particles. In 
the phases discussed in this article, the electric modes are 
spinless. 

^® Note that g and g~^ needn't be in the same conjugacy 
class. 

We have to be specific when we say that different defect 
condensates give different answers. For example, condens- 
ing I [123] > or I [124] > will give isomorphic answers, which 
is why we don't both put them up in our table. However, 
the TrS are different isomorphic subalgebras of A. 
One reads this off the magnetic part of the Hopf symmetry. 
If it is of the form F(G/H), where H is not a normal 
subgroup of G, then G/H is not a group, and one can prove 
that the algebra is not Hopf in this case. If the magnetic 
part is of the form F{Hm), where Hm is a group, then the 
residual symmetry algebra is a Hopf algebra. 



